arXiv: 1509.00037v2 [astro-ph.EP] 13 Oct 2015 


Spitzer IRAC Photometry for Time Series in Crowded Fields 


S. Calchi Novati^’^’^’“, A. Gould"^, J. C. Yee^’^, C. Beichman^, G. Bryden®, S. Garey^, 
M. Fausnaugh^, B. S. Gaudi*^, G. B. Henderson®’'^’'^, R. W. Pogge'^, Y. Shvartzvald®’'^, 

B. Wibking^ W. Zhu^ 

(Spitzer team) 
and 

A. Udalski^, R. Poleski^, M. Pawlak®, M. K. Szymanski®, J. Skowron®, P. Mroz®, 

S. Koziowski®, L. Wyrzykowski®, P. Pietrukowicz®, G. Pietrzynski®, I. Soszynski^, 

K. Ulaczyk^ 

(OGLE group) 

ABSTRACT 

We develop a new photometry algorithm that is optimized for Spitzer time 
series in crowded helds and that is particularly adapted to faint and/or heavily 
blended targets. We apply this to the 170 targets from the 2015 Spitzer mi- 
crolensing campaign and present the results of three variants of this algorithm 
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in an online catalog. We present detailed accounts of the application of this 
algorithm to two difficult cases, one very faint and the other very crowded. Sev¬ 
eral of Spitzer’s instrumental characteristics that drive the specihc features of 
this algorithm are shared by Kepler and WFIRST, implying that these features 
may prove to be a useful starting point for algorithms designed for microlensing 
campaigns by these other missions. 

Subject headings: gravitational lensing: micro - techniques: photometric 


Introduction 


The IRAC camera on the Spitzer observatory has been employed in a vast array of 
scientihc investigations, including some that were not imagined at the time it was launched 
in 2003, much less when it was conceived in the 1980s. Perhaps the most famous of 
these is the measurement of differential exo-planet transit depth as a function of wave- 
l ength, which requires exquisite photometric pr ecision in order to probe the exo-atmospheres 
f Charbonneau et al.ll2005l : iDeming et al.ll2005l) . These and other applications have stimu¬ 
lated the development of new photometric techniques and pipelines. 

In the last two years , Spitzer has taken on a completely new role as a “microlens p arallax 
satellite” flRefsdallll966l) . Although this application was sugges ted prior to laun ch fiGould 


1999h . and even carried out for one event fOGLE-2005-SMC-001. lDong et al.ll2007l) . the pho¬ 
tometric advances required to meet this role were not even basically understood until data 
were collected from the 2014 “pilot program”: 100 hours of Director’s Discretionary Time 
that were granted to assess the feasibility of Spitzer microlens parallaxes. These data revealed 
that while standard {Spitzer or other) photometric packages returned reasonably good pho 


tometry for most brig ht events (jUdalski et al.l l2015al: lYee et ahl l2015al: iGalchi Novati et ah 


20151: IZhu et al.ll2015al ). these packages usually failed, often catastrophically, for faint events. 


After 832 hours were awarded for the 2015 season (i.e., the majority of the 38 days 
for which Spitzer can observe Galactic bulge microlensing helds during microlensing season), 
Yee et ahl fj2015bl) developed detailed protocols to optimize efficiency to achieve the program’s 
primary aim of measuring the “Galactic distribution of planets”. Of course, one of the key 
parameters in these protocols was an es timate of the IRA G 3.6 fim threshold at which one 
could expect reliable photometry. I Yee et ahl fj2015bl) reviewed the results from various 

photometric packages as applied to the 2014 campaign to establish a threshold assuming no 
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further improvements, which they parameterized by a 3.6/rm flux indicatoi0 Leg = 15.5. 
They also analyzed the reasons for the shortcomings of existing packages when applied to 
this data set and reported in outline the ongoing work (now presented in the current paper) 
to resolve these shortcomings. 


One of the key points made by lYee et ahl fl2015bl) was that the selection of targets had 
to balance two critical indicators, encapsulated in a “quality factor”, which sums over all 
observed events, i, 


Q = }^SiP, 


( 1 ) 


where Si is the planet sensitivity and Pi is the probability that the observations would 
return a reliable microlens parallax (which of course requires reliable 5pzteer photometry). 
This implied that events with exceptional sensitivity (large Si) should be observed despite 
great risk (low Pi) that they would be too faint by the time Spitzer observed them, which 
could be from three to 10 days after they were recognized as sen sitive to planets based on 
their ground-based light curves (Figure 1 of lUdalski et al.ll2015al) . 


The main class of su ch ultra-sensitive events is high-magnihcation events flGriest fc Sahzadeh 


19981: iGould et al.ll2010h . which in most cases cannot be predicted until a few days before 
peak. Hence, in most cases, by the time that Spitzer observed such a high-magnihcation 
event it would already be well past peak and therefore very faint. Acting on the basis of 
Equation ([T]), however, the 5'pzfeer microlens team typically attempted to observe such high 
magni hcation events eve n when the expected 3.6 /xm hux was well below the threshold identi- 
hed bv lYee et ahl (I2015b|l . In fact, seven high-magnihcation events (with peak magnihcations 
A > 100) were observed, which greatly exceeded expectations. This also greatly increased 
the stakes of improving the photometric performance at faint hux levels. 

More generally, Fpxfeer microlens targets can be subject to extreme crowding, which can 
critically impact the photometry of these sources even when they are relatively bright. This 
can in turn adversely ahect a range of applications from m easuring microlens pa rallaxes to 
using the Spitzer light curves themselves to detect planets fiGould fc Hornel 1201311 


m 


The pro t ocol for the target selection for the 2015 Spitzer campaign is described in detail 


Yee et ahl (j2015bll . A key point is that Spitzer targets can be chosen ’’objectively” or 
’’subjectively”, according to which they fulhll, or not, a set of specihed criteria. This is 
relevant hrst to establish the cadence of the observations and second to frame the analysis 
of the planet sensitivity. 


^Leff = I — 0.93 A/ — 1.3 + 0.50(7 — Aj — 17.2), where 0 is a step function and Aj is the extinction. 
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In Section [2l we review the challenges to obtaining excellent time series photometry 
in crowded helds with Spitzer using existing packages, and we outline our algorithm for 
meeting these challenges. In Section [31 we discuss present applications of this algorithm 
to two examples of the difficult conditions encountered while pursuing the above scientific 
objectives. In Section 01 we describe an online catalog of photometry from our pipeline for 
all of the 170 events observed by Spitzer in the 2015 campaign. As we note there, it may 
often be possible to improve on this pipeline photometry for individual events of special 
scientific interest. The first author should therefore be contacted regarding applications of 
this photometry. Finally, in Section [5l we discuss some future implications of this work. 


2. Challenges to Crowded Field Photometry with Spitzer 


There are broadly two classes of routines applicable to time series photometry in crowded 
fields using electronic detectors: multi-object point spread function (PSF) photometry (MOPSF) 
and difference imaging analysis (DIA). 


MOPSF (e.g., DoPhot, ISchechter et ah! Il993l) was originally designed to disentangle 
stars in single-epoch images of crowded helds, such as globular clusters. Once an ensemble 
of point sources is identihed, they are ht either individually, or in groups to an ensemble 
of overlapping PSFs. In principle, this routine can be applied to a time series by simply 
repeating the operation on many frames and cross-identifying point sources. However, if 
the cadence is high (so that the stars can be approximated as not moving between epochs), 
then the stability of the results can be improved by enforcing hxed positions on the stars. 
This is particularly relevant for applications such as microlensing in which one star may 
substantially change brightness while the others remain roughly constant. In that case, if 
the star positions are not held hxed, they may drift under the inhuence of changing hux 
ratios of overlapping PSFs, including unresolved (i.e., unmodeled) blends. 


DIA 


e.g., 


Alard fc Luptonlll998l) works by a completely different principle. No attempt 


is made to identify stars. Rather, an essentially noise-free template is created by stacking 
good-seeing images. Each image in the series is then aligned photometrically and geometri¬ 
cally with the template, and the template is convolved with a kernel to make its PSF similar 
to that of the current image. Then the current image is subtracted from the convolved 
template. Stars that have varied then appear as PSF bumps or divots. Usually these are iso¬ 
lated and so can easily be centroided and then multiplied by the PSF to yield a photometric 
measurement. DIA has proved so far superior to MOPSF that the latter is rarely used for 
time series photometry in crowded helds. At the same time, however, it is obviously useless 
for photometry of constant stars since these simply disappear from the difference images. 
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Neither of these techniques is suitable for the analysis of faint variable targets in crowded 
Spitzer images. 

For the Spitze r microlensing campaign we use channel 1 (at 3.6/im) of the IRAC camera 
flFazio et ahl 120041 1 with mean pixel scale 1.221" and mean FWHM 1.66", with the image 
quality being limited primarily by the telescope optics. Given our goal of performing point- 
source photometry in a crowded held, it is important to note that the PSF is undersampled 
and tha t there is a signih cant variation in sensitivity within pixels (i.e., with the “pixel 


phase”, lingalls et ahl 1201211 . To address these issues the Spitzer Science Center developed 


an approach based on Point Response Functions (PRFs) specihc for IRAC (a tec hnique 
previously developed for the WFPC2 and NICMOS on the Hubble Space Telescope, iLauer 
19991: lAnderson &: Kindl2000h . These are meant to combine information about the PSF, 
the detector sampling, and the intrapixel sensit ivity variatioij^. The Spitz er Science Center 
has also developed a specihc package, MOPEX flMakovoz fc Marleaull200511 . to deal with the 
various aspects of the Spitzer instruments image analysis, and specihcally for the photometry 
of point sources for IRAC. 


We pause to note why the standard MOPSF and DIA packages are bound to fail, at 
least in the more difficult situations of low signal-to-noise ratio (S/N) and/or crowding. DIA 
does not work because it implicitly assumes a uniform pixel response, a condition that is 
not met for IRAC. MOPSF also does not work, but not for such a fundamental reason. The 
underlying logic of standard MOPSF is actually very similar to MOPEX, but it is not built 
to deal with the PRFs. In this sense MOPEX can be considered a specialized variant of 
MOPSF. 


For Spitzer microlensing observations, each “epoch” of a given target is composed by six 
30s exposures, each dithered by a few arcsec. The Spitzer pipeline (in particular, MOPEX) 
provides, in addition to these single frames, also a mosaic of the dithers of a given epoch. We 
remark that a PRF analysis only makes sense when applied to the single frame images. This 
information is blurred away in the mosaic frames, which come, on the other hand, with a 
higher S/N and look similar to optical CCD images. Standard photometry techniques, such 
as aperture photometry, MOPSF, and DIA can be applied to these images, and for bright 
and/or isolated stars the results can be passable to good. However, for faint stars in more 
crowded situations the results can be poor to catastrophic. 

Point-source photometry requires first, point-source extraction and second, flux deter¬ 
mination. In our analysis we take advantage of two facts, which then requires us to update 
and refine the MOPEX approach, accordingly. First, we possess quite accurate prior knowl- 


^http://irsa.ipac.caltech.edu/data/SPITZER/docs/dataanalysistools/tools/mopex/mopexusersguide/home/ 
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edge of essentially all significant 3.6 /xm sources on the frame. Second, we are dealing with a 
time series of images for which, as a working hypothesis, one can assume that only the mi- 
crolensed source varies, with all other neighboring stars assumed constant, including those 
that are heavily blended with the microlensed source. (In principle, it is possible that a 
neighboring star can be a variable. However, it is straightforward to identify such stars from 
the residuals, and to grant them an extra flux parameter for each epoch instead of imposing 
a single constant-flux parameter.) 

As foreshadowed, for point-source photometry, MOPEX works according to essentially 
the same principles as DoPhot, but with the specialized feature of employing PRFs. In 
essence, the photometry procedure is a £t based on a minimization where the signal, over 
the area of interest (a few pixels around each target), is modeled as a background term plus a 
linear sum of the PRFs of the different sources contributing to each pixel. There are therefore 
3 free parameters for each source: an amplitude term, which multiplies the corresponding 
PRF, and a pair of coordinates specifying the source centroid. As a technical point we note 
that the PRF is sampled at 5 x 5 = 25 pixel phase positions, i.e., five steps across the 
pixel in each direction. This is the basis for interpolating the PRF at each specific pixel 
phas^. MOPEX can perform this fitting scheme working on single frames individually or 
on multiple frames simultaneously, including frames belonging to different epochs. From our 
standpoint, the drawback of MOPEX is its rigidity: it does not allow us to take advantage 
of the two aforementioned facts specific to our problem. In particular, the only choice 
permitted by MOPEX is whether to leave the background term as a free parameter. The 
flux and positions of all the stars involved must always be fitted, which in our case means 
for each individual epoch (since the microlensed source flux changes with epoch). For this 
reason PRF-fitting MOPEX indeed works well with our data for sufficiently isolated and/or 
bright sources. However, especially for faint stars, it sometimes cannot locate the star at 
all, or if it can, incorrectly locates the centroid due to confusion of blended objects and 
so applies the wrong PRF. Or, what is potentially worse, as the microlensed source gets 
brighter and fainter, MOPEX shifts its estimate of the centroid, leading to systematic errors 
that vary with changing brightness. For the 2014 campaign, MOPEX worked reasonably 
well for about half of the sample. It worked poorly for about one quarter of the sample and 
failed catastrophically for the remaining quarter. MOPEX fared much worse on the 2015 
sample because the microlensed sources were overall substantially fainter. 


^In the current analysis we make use of the publicly available cold-mission PRFs. These are how- 
ever expected to have changed since 2009 when Spitzer entered its warm-mission operational phase 
( Storrie-Lombardi fc Dodd 2010[) . We have tested our algorithm using a preliminary version of the warm- 
mission PRFs (J. Ingalls, private communication) and we do not find substantial differences in the output 
photometry. 
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From an algorithmic standpoint, the path to resolving the difficnlties faced by MOPEX is 
straightforward: apply the same improvements made to DoPhot and other MOPSF packages 
when they gradnated from single-epoch to time-series applications. First, one shonld hold 
the brightness of all stars hxed from epoch to epoch (after establishing that they are not in 
fact variable). Second, one shonld (usnally) hold fixed the angnlar offsets between all stars 
from one epoch to the next. 

Implementation of the first (photometric) constraint is straightforward. If there are 
N epochs, n non-varying stars in the neighborhood of the lens, and m > 1 varying stars 
(inclnding the lens), then there are a total of n + mN + 1 photometric parameters (inclnding 
one for the backgronnd). 

However, it is not obvions how to implement the second (astrometric) constraint, or 
even whether to do so. At the ontset, prior to and independently from the photometric 
analysis, we bnild a relative astrometric solntion to link the OGLE (or in some cases, MOA) 
catalog coordinate^ to the Spitzer pixel coordinates, and this for each one of the QN IRAC 
frames. In one case (Method 2, below) we additionally need a relative astrometric solntion 
between an arbitrarily chosen IRAC “reference” frame and the other IRAC frames. 

The great majority of 3.6 fim sonrces in most microlensing helds are known with good 
astrometric precision from OGLE catalogs, which have mnch better resolntion than Spitzer. 
A first approach is simply to adopt these positions. In this case, we determine the positions 
of each star in each Spitzer frame, inclnding the microlensing sonrce, from a local ensemble 
of reference-star triangles. 

Second, one can nse the OGLE positions as inpnts to derive the IRAC positions from a 
simnltaneons fit to all frames. Then there are (6gA^-|-2(n-|-m)) parameters, i.e., some nnmber 
q > 6 frame-transformation parameters for each of the 6N images, which are applied prior to 
the analysis of each image, and the two-dimensional offsets of each star from a fidncial origin, 
which is held hxed for all 6N images. More specihcally, to carry ont the relative astrometry 
we make use of the IRAF package geomap with a polynomial non-linear transformation of 
order 5 to 8. The typical astrometric uncertainty that we achieve is about 0.1 pixel. 

Finally, one can use the OGLE positions as inputs, but derive star positions on each 
frame separately. From the standpoint of astrometry, this is similar to what MOPEX does 
on a single frame (although it still differs from MOPEX in that it imposes a photometric 


All Spitzer targets are selected from ongoing microlensing events that have been publicly alerted by either 
the Optical Gravitational Lens Experiment (OGLE) or Microlensing Observations in Astrophysics (MOA) 
collaborations. Each OGLE alert is linked to a detailed photometric and astrometric catalog of neighboring 
stars. 



constraint). In this case, there are 6N x 2{n + m) astrometric parameters. 
In brief, 

Method 1: Externally determined PRF centroid. 

Method 2: Single PRF centroid from simnltaneous £t to all epochs. 
Method 3: Individually htted PRF centroid for each frame. 


We hnd that in the majority of cases. Method 2 works best, but Methods 1 and 3 do 
sometimes work better. 

Specifically, if the microlensed source is extremely faint, there is almost no information 
in any of the images about the location of this source, so the OGLE information of this 
source position (which is very precise because it comes from astrometry of the difference 
image between when it is magnified and at baseline) is far superior to what can be obtained 
by fitting the IRAC images. Hence, Method 1 works best for these sources. 

If the microlensed source has considerable signal even in just a subset of images, then 
PRF hts to its position are usually superior to those derived from the transformation from 
the OGLE frame. This is an empirical fact. Rs cause is not completely understood but is 
most likely due to a combination of errors in the frame transformations and the different 
distributions of I band and 3.6/xm light in these crowded fields. Hence, either Method 2 or 
Method 3 is better. 

If the microlensed source is sometimes relatively bright and sometimes quite faint, then 
Method 2 works well but Method 3 fails because the lens position is “lost” on the fainter 
images. By the nature of the Spitzer microlensing campaign, this is true in the majority of 
cases. However, if the microlensed source is always bright, then Method 3 can work best, 
probably because it avoids the above-mentioned frame-transformation errors. 

In conclusion, our approach to the photometry for the microlensing campaign is a variant 
of MOPEX, a PRF-based minimization modeling of the signal^, where we have introduced 
the freedom to keep some parameters fixed from epoch to epoch, while other parameters are 
allowed to vary. Specifically, we always keep all other star fluxes fixed, except for those deter¬ 
mined to be variables. We usually keep all of the vector angular separations fixed (usually at 
values determined from a joint fit to the images, but sometimes by direct input from OGLE 


®In its present version for the fit we make use 
package https://cernlib.web.cern.ch/cernlib/ 


of MINUIT ( .Tames fc Roos 197fih within the CERNLIB 
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astrometry). For a few events that remain bright thronghont the Spitzer observations, we 
£t for the positions in each frame separately. 

In a nntshell, we take advantage of the PRF-htting procednre and at the same time of 
the specihc features of microlensing and of available data. 


3. Two Example Applications 

We illustrate our algorithm with two examples, Spitzer targets OGLE-2015-BLG-1395 
and OGLE-2015-BLG-0029. These are chosen because they present conditions that are 
challenging in substantially different ways. 


Applications of our algorithm to t wo other events ha ve already been published, the 


planetary event OGLE-2015-BLG-0966 ( 

Street et ah 

2015). and the binary event with a 

massive-remnant candidate OGLE-2015-BLG-1285 (£ 

hvartzvald et ah 

2015 

). We note that 


the case of OGLE-2015-BLG-0966 is particularly interesting because the early light curve 
yielded satifactory photometry with Method 1, but not Method 2. However, later, when the 
event became brighter. Method 2 became much better (although still not quite as good as 
Method 1) because it was able to apply its “knowledge” of the microlensed source position 
from the brighter images to perform precise photometry on the fainter ones. 

3.1. OGLE-2015-BLG-1395: High-Mag From Earth, Ultra-Faint Prom Spitzer 


OGLE-2015-BLG-1395 was discovered by the Optical Gravitational Lens Experiment 
(OGLE) based on observations with the 1.4 deg^ camera on its 1.3m Warsaw Telescope at the 
Las Gampanas Obs ervatory in Ghile using its Ea.rlv W arning System (EWS) real-time event 
detection software f Udalski et ah 1994 : Udalski 2003 ) on 19 June UT 17:56. The Spitzer 
team contacted OGLE 21 hours later to get an early update on the next night of data and 
based on this immediately (20 June UT 14:30) alerted the r nicrolensing community that this 
was a high-mag event. For results of this monitoring, see iGole et al.l (120151) . O n 22 June, 


the eve nt met the “objective criteria” for selection according to the protocols of lYee et ah 


fj2015b[) and was scheduled for observations at a cadence of 2/day for the hrst week and 
1/day thereafter. The hrst S'pzUer observations were three days later, HJD 7199.50, i.e., 
when the ground-based event was 5.3 days past peak and the event had already fallen to 
I = 18.6. In fact, when the Spitzer data are aligned to the OGLE scale (see below) the hrst 
point is at / = 19-0. cor responding to Lefr = 16.2, roughly 0.7 mag fainter than threshold 
set bv lYee et al.l fj2015bll . Hence, it was uncertain whether the parallax of this event could 


be recovered. 
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Figure [T] shows four reductions of the Spitzer data. MOPEX clearly fails. The reason 
for this failure is that it “finds” the source at ~ 3-4" from its true position. Regarding the 
three variants of the new reductions, Method 1 is satisfactory, Method 2 works best, and 
Method 3 fails completely. 


For many years it was believed that satellite parallax measurements required that the 
satellite observe the event over peak, or at l east close enough to peak to dete rmine the 
satellite-based impact parameter, uo.sat (e.g., iGouldl Il995l: iGaudi fc Gouldl 119971) . This is 
because satellite-based parallaxes tte are basically determined from 


tte = ^ 13 )] 3 \t = ; A/3 = ±no,e - ±Mo,sat, (2) 

Ui_ ie 

where the subscripts indicate parameters as measured from Earth and the satellite. Here, 
(fo,Wo,tE) are the standard point-lens microlensing parameters: time of maximum, impact 
parameter, and Einstein timescale, while is the Earth-satellite separation projected on 
the sky. 


However, lYee et ahl (j2015bl) argued that if the Spitzer source flux could be indepen¬ 
dently measured (rather than fit from the Spitzer light curve as was previously believed 
necessary), then tg.sat and Mo,sat could be determined even if the Spitzer data covered only 
a fragment of a post-peak light curve by making use of this independent determination 
of fs together with the value of ^e measured from Earth. From Figure [H it is clear that 
OGLE-2015-BLG-1395 provides an excellent way to test this idea. 


To implement it, we first derive an estimate of the (/ogle — Tspitzer)® color and then 
incorporate this directly into a joint fit to OG1jE+S pitzer data. The general method for 
estimating (/ — L)s is to measure the source color in some ground-based bands X — I, then 
to determine an XIL {X — I vs I — L) color-color relation based on cross-identified field 
stars in ground and Spitzer images, and hnally to apply this relation to the {X — L)s color 
to derive (/ — L)s. For example, OGLE often has sufficient H-band data during the event to 
measure V — I from the ground-based light curve. Alternatively, if the ht to ground-based 
data shows that the source is unblended, then the baseline [V — I) color can be used as 
the source color. However, in the present case, the event peaked too rapidly to determine 
{y — I)s from the normal cadence of OGLE H-band data. 

We therefore use the simultaneous I and H data taken with the dual-channel ANDIGAM 
camera on the 1.3m SMARTS telescope at GTIO, which acquired substantial data on almost 
all 2015 Spitzer microlensing events, primarily for exactly this purpose. Figure [2] shows a 
model-independent regression of these data, which yield (/ — //)s,ctio = —0.624 ± 0.007. We 
then form an IHL color-color diagram (Figure [3]). We carry out the color-color regression on 
stars —0.5 < (/ — //)ctio < +0.2 in order to restrict the sample to warm (T > 4300 K) stars 
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that are in the bulge (and so behind the same column of dust), and hnd Jctio — -hspitzer = 
1.210(1 — i/)ctio — 3.196 and hence (/cuo — -f^spitzer)* = —3.95. By regression, we obtain 
-fctio - -fogie = 0.73 and so hnally (Jogie - Tv^pitzer)* = -4.68. Here L = 25 - 2.5 log(Fspit2er) 
where Fgpitzer is the instrumental flux shown in Figure [T] (three upper panels). We estimate 
an uncertainty of ±0.1 on this transformation. We report below on how the final results 
change with the adopted error bar. In both cases we exclude recursively the outliers, also 
indicated in the figure, from the regression analysis. 

We apply this color constraint to a joint fit to OGLE and Spitzer data (Fig. 0]), wherein 
the statistical error bars are renormalized so that x ^/dof = 1 for each observatory. This is the 
standard procedure to account any possible additional scatter due to unmodeled systematics 
(this issue is further addressed in Sect. 0]). Specihcally, the size of the error bars is about 
0.08 and 0.12 mag at the bright and the faint end, respectively. Actually, there are four 
possible fits, parameterized by ±mo,® and ±Mo,sat as indicated by Equation ([2]). From the 
present standpoint of illustrating how Spitzer photometry is carried out and how external 
constraints are applied to the light curve, (to, wo)sat are the primary quantities of interest for 
these solutions. We find. 


(fo,Mo)++ = (7197.50, ±0.280); (to, Mo)+-= (7197.60, ±0.282) (3) 

(to, Mo)— = (7197.48,-0.280); (to, mo)-+ = (7197.60,-0.282), (4) 


where the hrst and second subscripts refer to the signs of Uq for Earth and Spitzer, respec¬ 
tively. That is, the naive idea that what is constrained by the Spitzer light curve is {to, |Mo|)sat 
is confirmed, even though the light curve does not cover the peak. However, this is true only 
because of the color constraint. If the uncertainty on the color constraint is raised from 0.1 
to 0.2, the above results barely change. However, if it is removed, then neither to nor uo is 
meaningfully constrained. 


These results confirm the conjecture of lYee et ahl (j2015bl) on the utility of obtaining 
post-peak light curves from Spitzer. In particular, they imply that planet-sensitive high- 
magnihcation events should be “subjectively selected” even when it is clear that Spitzer 
observations cannot begin until several days after peak. 


3.2. OGLE-2015-BLG-0029: Bright Source in a Crowded Field 

OGLE-2015-BLG-0029 was discovered by OGLE on 12 February 2015. On 10 May it 
was “subjectively selected” for Spitzer observations with “objective” 1/day cadence. This 
was long before Spitzer observations began, which is critical for establishing the planet 
sensitivity of events that do not meet objective criteria. According to the protocols of 
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Yee et al.l fl2015bl) . planets and planet sensitivity for snbjectively selected events begin with 


the time of selection, not the start of Spitzer observations. As it tnrned o nt, on 1 Jnne, the 
event met the objective criteria for rising events as laid ont in lYee et al.l fl2015bl) . Since the 
cadence for Spitzer observations of objectively selected events was the same as that specihed 
on 10 May, Spitzer observations of OGLE-2015-BLG-0029 were schednled objectively. These 
observations began on 11 Jnne, when the event hrst became visible by Spitzer^ and continned 
nntil the end of the campaign on 19 Jnly. As specihed in I Yee et ahl (j2015bl l. a cadence of 
“1/day” determines the relative freqnency with which an event should be observed, and 
the actual cadence is set by this and the available time. Hence, this event received extra 
observations at the beginning and end of the campaigns, as competing events were eliminated 
due to Spitzer's Sun-exclusion angle. Altogether, it was observed at 52 epochs. 


In contrast to OGLE-2015-BLG-1395, OGLE-2015-BLG-0029 is a very bright event with 
I ~ 14.5-15. However, similar to OGLE-2015-BLG-1395, the MOPEX pipeline basically 
fails to measure the light curve for this event, albeit not catastrophically, as was the case 
for OGLE-2015-BLG-1395. See Figure [5l The reason that MOPEX fails is crowding, in 
particular the presence of an / = 17 neighbor at 0.97". This neighbor, which is ~ 5 times 
fainter than the microlensed source at 3.6/xm at its faintest (hrst) Spitzer epoch, is clearly 
resolved in OGLE data but unresolved in individual IRAG images. What makes this case 
especially instructive is that the source lies almost exactly at the hrst Airy null (1.06"), so that 
it is resolvable in principle, but within one IRAG pixel (1.2"). If MOPEX is “informed” about 
the presence of this source, then it assigns the source wildly varying positions and huxes, 
which impacts the microlensing light curve catastrophically. If MOPEX is not ’’informed” 
(and so treats the microlensed source and the nearby blend as having a common position), 
it does better but still has more than an order of magnitude greater scatter than our new 
algorithm, particularly in its Method 2 variant. 


4. Description of Catalog 


In Table [T] we list the 170 events monitored in 2015. For each, we report the event 
name, the coordinates, the hrst and last day of observation, and the number of observed 
epochs. The events w ere chosen based on the micr olensing alerts provided by the OGLE 
fjUdalski et al.l l2015bl) and MOA flBond et ahl 120041) collaborations. The current analysis 
is based on the preliminary reduced data made available by the Spitzer Science Genter 
almost in real time (on average, 2-3 days after the observations). The hnal reduction of 
the data is now publicly available at the NASA/IPAG Infrared Science Database (IRSA, 
http://irsa.ipac.caltech.edu/frontpage/). Besides the full table, in the online material we 
provide the light curves for the full sample of observed events, similar to Figures [T] and 
E] (excluding of the MOPEX reductions). Finally, an updated version of our photometry 
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solutions is kept at a publicly available websitqj. In particular we provide tables with the 
data for all three reductions discussed in the text. The flux is expressed in arbitrary units. In 
some cases, none of the reductions look satisfactory, which may just reflect that the absolute 
level of source flux variations is very low. The reported error is the purely statistical error 
from the fitting procedure based on the IRAC uncertainty images (and so underestimates the 
true uncertainties). In particular, we attribute any additional scatter to systematics due to 
unmodeled physics (a completely generic feature for all photometry algorithms in crowded 
fields). These must be considered case per case during the modeling of the events. We 
note that a single parameter (the “background”) must account for the combined effects of 
unresolved stars and incompletely modeled wings of distant bright stars, in addition to the 
true background. It is therefore expected that faint stars will sometimes take on negative 
flux values. This has no impact on microlensing analyses, which £t only the flux variations 
and similarly absorb the constant-flux term into a nuisance parameter “/f,”. 


We emphasize that for events of special scientific interest, further work on reductions 
may improve the light curve. In particular, for events in which the external knowledge of 
the source position is important, such as OGLE-2015-BLG-1395, this can often be improved 
by measuring this position from OGLE difference images. An additional key aspect of 
the analysis that may be improved on a case by case basis is the choice of blend stars 
to be fitted simultaneously with the microlensed source. An additional aspect which may 
become critical, which is not taken into account in the standard pipeline and that requires 
further specific analyses, is saturation and nonlinearities in the response of the detector 
when approaching the pixel full well. An example is that of the Spitzer light curve for the 
microlensing event OGLE-2015-BLG-0763, which is presented in the online catalog without 
the special treatment required to deal with saturation. Because this event is of special 
interest, it w a.s reanalyzed mar king use of specific procedures to deal with saturation and 
non-linearity (jZhu et al.ll2015bll . 


Discussion 


The principal characteristics of the data set th at drove the desig n of the algorithms 
described here will also be present in the Kepler K2 flHowell et al.l 120141 ) microlensing data. 
These include. 


1) Undersampled PSF, 

2) Non-uniform pixel response. 


® htt p: / / WWW. astronomy, ohio-state .edu/Spitzer2015/ 
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3) Many stars overlapping the microlensed sonrce, 

4) Most (or all) of the other stars non-variable, 

5) External astrometric catalog at higher resolution, 

6) Dithered exposures, 

7) Many epochs. 


Hence, it may be prohtable to employ similar algorithms on the K2 data set. While the hrst 
microlensing data will not be available at least until June 2016, there are already K2 test 
data from other crowded helds on which these algorithms could be developed. 


All of the above ch aracteristics, with th e exception of (5) and the partial exception of 
(3) apply to WFIRST fISpergel et al.ll2013alJbl ) as well. In some respects, the challenges of 
WFIRST may seem less pressing. However, if having available a source catalog at higher 
resolution proves important for WFIRST (as it so far appears to be for Spitzer), then now 
is the time to remedy its absence, before the demise of the Hubble Space Telescope. 


We thank J. Ingalls for a useful discussion about PRF htting with IRAC data. Work 
by SCN, AG, SC, JCY, and WZ was supported by JPL grant 1500811. Work by YS was 
supported by an appointment to the NASA Postdoctoral Program at the Jet Propulsion Lab¬ 
oratory, administered by Oak Ridge Associated Universities through a contract with NASA. 
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(Caltech)/Jet Propulsion Laboratory (JPL) funded by NASA through the Sagan Fellowship 
Program executed by the NASA Exoplanet Science Institute. This work is based (in part) 
on observations made with the Spitzer Space Telescope, which is operated by the Jet Propul¬ 
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Fig. 1.— Spitzer photometry of OGLE-2015-BLG-1395 using four different packages: 
MOPEX, and the three variants of the algorithm presented in this paper; Method 1: ex¬ 
ternally hxed centroid, Method 2; centroid jointly fitted to entire light curve. Method 3: 
centroid htted separately for each epoch. Both MOPEX and Method 3 fail completely and 
for the same reason: very faint target. 
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F| (Instrumental) 


Fig. 2.— Model-independent regression of iJ-band on /-band instrumental fluxes from 
CTIO-SMARTS observations of OGLE-2015-BLG-1395. These yield a model-independent 
instrumental source color {I — H)s,ctio = —2.5 log(F//Fff) ~ —2.5 log(4.00/2.25) = —0.624 ± 
0.007. The circles indicate the points used for the color determination whereas recursively 
excluded points are indicated by crosses. 




(I-L) (instrumental) 
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(l-H) (instrumental) 


Fig. 3.— Color-color diagram Jctio — -^spitzer vs. (/ — Fr)ctio based on field stars near OGLE- 
2015-BLG-1395. Only stars —0.5 < (J — ir)ctio < 0.2 are included in the fit to ensure that 
it is dominated by bulge stars that are reasonably warm T > 4250. The open points have 
been recursively excluded. The blue line is the source color measured in Figure [2l From this 
we derive an instrumental color (Jctio — J^spitzer)® = —3.95 ± 0.10 . 
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>190 7200 7210 7220 

HJD - 2450000 


Fig. 4.— Joint point-lens parallax fit of OGLE and Spitzer observations of OGLE-2015- 
BLG-1395, nsing Method 2 rednctions for Spitzer. The OGLE data (black) are shown at 
their observed calibrated J-band magnitnde. The Spitzer data (red) are scaled so that 
their “/ magnitnde” corresponds to the same magnihcation as an OGLE point at this same 
magnitnde. The ht is constrained by the (/ogle—-hspitzer)s = —4.68±0.10 sonrce color, which is 
taken from Fignre[3] (after acconnting for the instrnmental color relation Jctio — Ggie = 0.68). 
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Fig. 5.— Four reductions of OGLE-2015-BLG-0029 Spitzer data. From top to bottom they 
are Methods 1, 2, 3, and the standard MOPEX pipeline. Although the microlensed source 
is bright, it is blended with a ~ 5 times fainter star at 1", which is similar to the first Airy 
null (1.06") but smaller than an IRAG pixel (1.2"). Method 2 copes with this situation 
best. MOPEX becomes completely confused if it is “informed” of the existence of this blend 
(not shown) and does better if it treats the blend as coincident with the microlensed source 
(displayed), but the results are still quite poor when compared to Method 2. 






































































Table 1. Observed Events 


Event 

RA (J2000) 

(degrees) 

DEC (J2000) 

(degrees) 

first 

(HJD) 

Spitzer 

last epochs 

(HJD) 

MOA-BLG-2015-020 

268.219917 

-32.485861 

7183.95 

7221.76 

61 

MOA-BLG-2015-079 

268.085334 

-31.295665 

7183.66 

7221.76 

112 

MOA-BLG-2015-117 

268.792652 

-35.080247 

7184.10 

7212.60 

45 

MOA-BLG-2015-220 

271.832485 

-25.467555 

7187.12 

7222.96 

48 

MOA-BLG-2015-237 

269.467804 

-31.472354 

7184.97 

7223.06 

54 

MOA-BLG-2015-240 

270.585946 

-31.996441 

7185.80 

7212.75 

37 

MOA-BLG-2015-245 

269.321669 

-27.857411 

7184.99 

7222.60 

78 

MOA-BLG-2015-267 

269.214339 

-30.775359 

7199.49 

7222.65 

39 

MOA-BLG-2015-282 

272.572708 

-29.265556 

7199.71 

7222.94 

29 

MOA-BLG-2015-287 

270.125760 

-33.017623 

7199.35 

7222.89 

36 

OGLE-BLG-2014-0298 

268.565625 

-31.018083 

7185.73 

7221.84 

25 

OGLE-BLG-2014-0613 

268.490300 

-28.572670 

7183.98 

7221.81 

158 

OGLE-BLG-2015-0011 

269.217833 

-29.283250 

7184.96 

7222.59 

53 

OGLE-BLG-2015-0022 

268.090792 

-29.785028 

7183.67 

7221.05 

60 

OGLE-BLG-2015-0029 

269.944167 

-28.644944 

7185.31 

7222.89 

52 

OGLE-BLG-2015-0034 

270.580333 

-27.516083 

7186.01 

7222.92 

62 

OGLE-BLG-2015-0060 

269.993125 

-27.780944 

7185.78 

7222.90 

50 

OGLE-BLG-2015-0081 

268.653000 

-28.996278 

7184.10 

7221.82 

57 

OGLE-BLG-2015-0145 

270.177833 

-35.154056 

7185.29 

7222.88 

51 

OGLE-BLG-2015-0149 

270.288208 

-32.557694 

7185.80 

7222.93 

49 

OGLE-BLG-2015-0195 

266.790792 

-33.417278 

7182.39 

7220.58 

70 

OGLE-BLG-2015-0196 

266.492958 

-32.956778 

7182.39 

7220.56 

70 

OGLE-BLG-2015-0211 

262.359083 

-30.981750 

7180.20 

7216.80 

109 

OGLE-BLG-2015-0238 

265.708250 

-36.818500 

7182.37 

7219.17 

75 

OGLE-BLG-2015-0244 

268.240625 

-28.980889 

7183.97 

7221.75 

59 

OGLE-BLG-2015-0254 

271.553292 

-26.850194 

7186.87 

7222.96 

49 

OGLE-BLG-2015-0326 

274.074625 

-26.338139 

7199.73 

7223.03 

29 

OGLE-BLG-2015-0350 

268.248583 

-31.820278 

7183.95 

7221.76 

61 

OGLE-BLG-2015-0379 

269.104292 

-29.574056 

7184.61 

7222.58 

54 

OGLE-BLG-2015-0388 

268.468917 

-28.534028 

7183.98 

7221.81 

66 

OGLE-BLG-2015-0401 

264.081167 

-27.774139 

7180.19 

7213.03 

94 

OGLE-BLG-2015-0444 

262.917375 

-30.673111 

7180.20 

7216.79 

84 

OGLE-BLG-2015-0448 

272.559917 

-31.752611 

7187.54 

7223.04 

211 

OGLE-BLG-2015-0461 

270.043208 

-28.156944 

7185.79 

7222.90 

58 

OGLE-BLG-2015-0477 

273.200500 

-25.014028 

7188.88 

7222.97 

45 

OGLE-BLG-2015-0479 

265.918958 

-35.509278 

7182.38 

7219.17 

67 

OGLE-BLG-2015-0529 

270.264667 

-29.922917 

7185.79 

7222.92 

51 

OGLE-BLG-2015-0561 

267.224542 

-29.674444 

7183.09 

7220.61 

68 

OGLE-BLG-2015-0565 

269.153708 

-29.128056 

7184.62 

7222.59 

53 

OGLE-BLG-2015-0572 

272.847750 

-25.473861 

7188.88 

7222.97 

44 

OGLE-BLG-2015-0607 

274.049625 

-27.451306 

7189.13 

7223.03 

54 

OGLE-BLG-2015-0642 

263.539458 

-24.882528 

7192.99 

7216.78 

25 

OGLE-BLG-2015-0692 

268.079708 

-28.133917 

7185.70 

7221.07 

73 
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Table 1—Continued 


Event 

RA (J2000) 

(degrees) 

DEC (J2000) 

(degrees) 

first 

(HJD) 

Spitzer 

last 

(HJD) 

epochs 

OGLE-BLG-2015-0703 

269.048583 

-29.853389 

7184.47 

7222.17 

65 

OGLE-BLG-2015-0709 

266.113292 

-34.442972 

7182.38 

7219.16 

74 

OGLE-BLG-2015-0713 

272.660042 

-28.412083 

7213.65 

7223.01 

28 

OGLE-BLG-2015-0749 

264.955708 

-27.511722 

7180.78 

7206.48 

59 

OGLE-BLG-2015-0761 

272.084667 

-25.443056 

7187.52 

7222.97 

46 

OGLE-BLG-2015-0763 

263.097542 

-29.302694 

7180.19 

7216.78 

142 

OGLE-BLG-2015-0769 

265.670083 

-27.368833 

7182.40 

7219.15 

59 

OGLE-BLG-2015-0772 

269.421583 

-28.143139 

7184.99 

7222.60 

54 

OGLE-BLG-2015-0798 

271.680583 

-27.590167 

7186.86 

7222.95 

47 

OGLE-BLG-2015-0802 

266.514625 

-36.497917 

7182.38 

7220.57 

76 

OGLE-BLG-2015-0808 

267.184833 

-33.102278 

7199.31 

7220.59 

25 

OGLE-BLG-2015-0828 

267.132917 

-34.184778 

7183.08 

7220.58 

72 

OGLE-BLG-2015-0843 

271.197042 

-27.176472 

7186.87 

7222.95 

61 

OGLE-BLG-2015-0845 

271.088708 

-31.580556 

7186.24 

7222.93 

148 

OGLE-BLG-2015-0877 

273.676958 

-27.970917 

7188.92 

7212.92 

29 

OGLE-BLG-2015-0912 

265.957583 

-26.519944 

7182.41 

7219.15 

57 

OGLE-BLG-2015-0914 

269.311750 

-29.760139 

7184.96 

7222.58 

53 

OGLE-BLG-2015-0925 

270.259292 

-26.384889 

7199.67 

7222.91 

29 

OGLE-BLG-2015-0926 

270.883792 

-26.135806 

7199.67 

7222.91 

29 

OGLE-BLG-2015-0930 

264.393917 

-25.512667 

7180.21 

7218.07 

111 

OGLE-BLG-2015-0941 

272.348083 

-25.762056 

7187.53 

7222.98 

47 

OGLE-BLG-2015-0944 

265.258000 

-27.853806 

7181.09 

7198.12 

42 

OGLE-BLG-2015-0955 

273.709833 

-27.006833 

7213.88 

7223.01 

14 

OGLE-BLG-2015-0958 

267.754958 

-28.416861 

7183.32 

7221.06 

64 

OGLE-BLG-2015-0961 

268.146375 

-30.080361 

7185.69 

7221.75 

131 

OGLE-BLG-2015-0965 

269.291083 

-28.961056 

7184.99 

7222.59 

52 

OGLE-BLG-2015-0966 

268.754250 

-29.047111 

7192.64 

7222.14 

119 

OGLE-BLG-2015-0968 

268.985833 

-28.854083 

7184.47 

7222.19 

56 

OGLE-BLG-2015-0977 

272.061375 

-26.159472 

7187.53 

7222.96 

47 

OGLE-BLG-2015-0984 

263.875167 

-26.553778 

7185.84 

7217.51 

39 

OGLE-BLG-2015-0987 

268.193708 

-28.361222 

7183.97 

7221.07 

65 

OGLE-BLG-2015-1020 

263.050417 

-29.781722 

7180.19 

7216.78 

108 

OGLE-BLG-2015-1076 

270.744792 

-32.853694 

7192.81 

7222.94 

135 

OGLE-BLG-2015-1077 

274.478083 

-22.992611 

7192.93 

7223.03 

64 

OGLE-BLG-2015-1084 

271.146417 

-26.168028 

7186.87 

7212.80 

35 

OGLE-BLG-2015-1095 

267.693875 

-32.708167 

7183.24 

7221.04 

65 

OGLE-BLG-2015-1096 

267.872333 

-32.153694 

7183.40 

7221.04 

62 

OGLE-BLG-2015-1100 

263.669042 

-29.885194 

7180.20 

7217.52 

84 

OGLE-BLG-2015-1109 

270.213083 

-30.645139 

7185.80 

7191.80 

13 

OGLE-BLG-2015-1112 

273.535250 

-25.482778 

7188.79 

7222.98 

51 

OGLE-BLG-2015-1113 

264.277208 

-27.299083 

7180.18 

7218.10 

88 

OGLE-BLG-2015-1123 

273.648042 

-26.507389 

7188.88 

7223.01 

70 

OGLE-BLG-2015-1124 

273.253750 

-27.532861 

7188.89 

7206.39 

23 
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Table 1—Continued 


Event 

RA (J2000) 

(degrees) 

DEC (J2000) 

(degrees) 

first 

(HJD) 

Spitzer 

last epochs 

(HJD) 

OGLE-BLG-2015-1129 

266.413125 

-27.160028 

7182.41 

7185.33 

18 

OGLE-BLG-2015-1136 

268.557583 

-28.708361 

7185.71 

7198.40 

24 

OGLE-BLG-2015-1145 

269.451042 

-29.080444 

7185.63 

7222.84 

79 

OGLE-BLG-2015-1148 

270.259625 

-28.686111 

7185.79 

7198.57 

22 

OGLE-BLG-2015-1150 

273.049208 

-25.661694 

7199.72 

7223.05 

121 

OGLE-BLG-2015-1153 

264.317792 

-29.426389 

7193.01 

7218.08 

50 

OGLE-BLG-2015-1161 

267.618375 

-29.979222 

7192.59 

7220.61 

34 

OGLE-BLG-2015-1167 

269.024500 

-28.626361 

7184.46 

7222.20 

56 

OGLE-BLG-2015-1172 

268.545417 

-31.103194 

7185.73 

7198.41 

24 

OGLE-BLG-2015-1184 

268.514458 

-30.311250 

7183.96 

7185.41 

10 

OGLE-BLG-2015-1187 

268.657792 

-31.681222 

7199.45 

7222.16 

27 

OGLE-BLG-2015-1188 

268.877125 

-31.281667 

7184.34 

7222.17 

56 

OGLE-BLG-2015-1189 

270.735000 

-30.373056 

7192.80 

7222.93 

87 

OGLE-BLG-2015-1194 

261.049625 

-28.947472 

7207.44 

7215.02 

9 

OGLE-BLG-2015-1196 

266.466667 

-24.287278 

7182.41 

7198.06 

37 

OGLE-BLG-2015-1197 

264.229917 

-29.782028 

7207.42 

7218.09 

31 

OGLE-BLG-2015-1204 

268.927083 

-29.384306 

7199.50 

7222.18 

27 

OGLE-BLG-2015-1207 

267.709667 

-23.445500 

7192.93 

7220.56 

31 

OGLE-BLG-2015-1212 

268.103292 

-29.181111 

7185.70 

7206.09 

92 

OGLE-BLG-2015-1221 

268.418792 

-32.279306 

7199.33 

7221.94 

101 

OGLE-BLG-2015-1223 

265.690333 

-27.528333 

7182.40 

7212.95 

53 

OGLE-BLG-2015-1227 

268.327417 

-29.652833 

7199.41 

7221.83 

26 

OGLE-BLG-2015-1232 

267.641125 

-32.739639 

7185.66 

7191.98 

16 

OGLE-BLG-2015-1234 

267.179917 

-24.520028 

7192.93 

7220.56 

31 

OGLE-BLG-2015-1236 

265.421375 

-27.809278 

7185.81 

7219.15 

42 

OGLE-BLG-2015-1241*" 

269.529708 

-27.932000 

7185.63 

7198.49 

35 

OGLE-BLG-2015-1256 

269.537083 

-28.769083 

7185.77 

7206.27 

31 

OGLE-BLG-2015-1262 

267.523792 

-30.712611 

7185.68 

7206.04 

33 

OGLE-BLG-2015-1263 

267.717208 

-32.380389 

7185.67 

7206.03 

32 

OGLE-BLG-2015-1264 

264.976250 

-34.584722 

7192.51 

7219.17 

31 

OGLE-BLG-2015-1268 

269.207375 

-21.899333 

7192.61 

7212.93 

48 

OGLE-BLG-2015-1281 

273.034625 

-27.167694 

7192.86 

7222.95 

36 

OGLE-BLG-2015-1285 

264.846875 

-27.820278 

7199.81 

7218.09 

39 

OGLE-BLG-2015-1288 

267.516625 

-33.368472 

7192.52 

7221.03 

35 

OGLE-BLG-2015-1289 

269.461583 

-28.217306 

7206.94 

7222.83 

22 

OGLE-BLG-2015-1295 

270.226042 

-27.302139 

7199.39 

7223.08 

127 

OGLE-BLG-2015-1297 

268.087542 

-29.683806 

7192.59 

7206.05 

17 

OGLE-BLG-2015-1303 

273.970333 

-30.373194 

7199.42 

7223.02 

68 

OGLE-BLG-2015-1309 

264.693250 

-27.479472 

7199.82 

7218.10 

20 

OGLE-BLG-2015-1315 

266.672583 

-33.458194 

7206.71 

7220.57 

17 

OGLE-BLG-2015-1319 

269.443333 

-32.472194 

7206.95 

7223.07 

22 

OGLE-BLG-2015-1325 

266.423625 

-35.957278 

7213.41 

7220.83 

18 

OGLE-BLG-2015-1328 

267.063208 

-32.687778 

7199.32 

7206.02 

8 
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Table 1—Continued 


Event 

RA (J2000) 

(degrees) 

DEC (J2000) 

(degrees) 

first 

(HJD) 

Spitzer 

last epochs 

(HJD) 

OGLE-BLG-2015-1339 

263.321917 

-30.271222 

7199.91 

7217.52 

18 

OGLE-BLG-2015-1341 

268.563667 

-28.350167 

7199.47 

7221.81 

25 

OGLE-BLG-2015-1344 

268.934542 

-27.940028 

7192.67 

7222.32 

57 

OGLE-BLG-2015-1346 

269.806958 

-28.936611 

7192.79 

7222.88 

37 

OGLE-BLG-2015-1348 

269.487708 

-30.509417 

7199.55 

7223.06 

29 

OGLE-BLG-2015-1350 

266.685583 

-36.592889 

7199.29 

7220.57 

73 

OGLE-BLG-2015-1352 

265.141875 

-27.338306 

7207.37 

7218.08 

13 

OGLE-BLG-2015-1366 

268.289250 

-21.944056 

7213.90 

7221.74 

10 

OGLE-BLG-2015-1368 

267.945083 

-29.575167 

7199.34 

7221.06 

77 

OGLE-BLG-2015-1370 

268.956667 

-29.124111 

7206.90 

7222.18 

19 

OGLE-BLG-2015-1371 

268.665875 

-29.093944 

7206.87 

7221.83 

18 

OGLE-BLG-2015-1374 

271.134375 

-26.887806 

7213.64 

7223.00 

28 

OGLE-BLG-2015-1378 

264.834792 

-22.976222 

7199.79 

7218.06 

20 

OGLE-BLG-2015-1383 

268.651292 

-28.702250 

7213.53 

7221.94 

22 

OGLE-BLG-2015-1392 

264.095125 

-27.298917 

7213.68 

7217.82 

11 

OGLE-BLG-2015-1395 

271.783917 

-26.303333 

7199.51 

7222.96 

36 

OGLE-BLG-2015-1400 

268.109625 

-28.866250 

7213.47 

7221.55 

20 

OGLE-BLG-2015-1403 

267.882875 

-29.340361 

7206.79 

7221.06 

17 

OGLE-BLG-2015-1412 

270.124875 

-27.435361 

7207.00 

7223.05 

35 

OGLE-BLG-2015-1416 

270.937500 

-26.848056 

7207.07 

7222.92 

21 

OGLE-BLG-2015-1420 

271.477500 

-28.557361 

7213.65 

7223.02 

28 

OGLE-BLG-2015-1424 

264.063292 

-29.642833 

7199.83 

7218.09 

20 

OGLE-BLG-2015-1425 

263.142958 

-30.716333 

7213.69 

7217.10 

9 

OGLE-BLG-2015-1435 

269.842167 

-27.219222 

7199.57 

7206.27 

8 

OGLE-BLG-2015-1440 

270.973167 

-27.707417 

7213.64 

7222.99 

28 

OGLE-BLG-2015-1447 

269.233417 

-31.941111 

7213.50 

7222.81 

27 

OGLE-BLG-2015-1448 

270.939333 

-27.852194 

7207.04 

7222.99 

35 

OGLE-BLG-2015-1450 

270.210000 

-28.671028 

7213.51 

7223.07 

56 

OGLE-BLG-2015-1457 

267.857167 

-30.436389 

7206.78 

7221.05 

17 

OGLE-BLG-2015-1467 

265.107500 

-23.801389 

7207.38 

7212.97 

7 

OGLE-BLG-2015-1469 

265.266583 

-24.848944 

7213.67 

7218.92 

13 

OGLE-BLG-2015-1470 

268.501833 

-31.592806 

7206.82 

7221.84 

18 

OGLE-BLG-2015-1481 

268.008083 

-30.985028 

7213.48 

7221.54 

20 

OGLE-BLG-2015-1482 

267.630542 

-30.888694 

7206.73 

7221.04 

17 

OGLE-BLG-2015-1485 

265.380250 

-24.214972 

7207.35 

7218.92 

20 

OGLE-BLG-2015-1492 

268.686000 

-29.707361 

7213.48 

7221.96 

22 

OGLE-BLG-2015-1500 

272.021583 

-27.445611 

7207.12 

7212.84 

7 

OGLE-BLG-2015-1521 

267.485792 

-21.411889 

7213.66 

7220.82 

17 

OGLE-BLG-2015-1530 

269.018958 

-28.240472 

7213.53 

7222.38 

24 

OGLE-BLG-2015-1533 

269.048083 

-29.476417 

7213.51 

7222.39 

24 

OGLE-BLG-2015-1553 

265.857042 

-23.377750 

7213.67 

7219.14 

14 


Note. — We report name, target coordinates and information about Spitzer obser¬ 
vations for the 170 microlensing events monitored during the 2015 Spitzer microlens- 







ing campaign. We note that 2 out of the 10 MOA events have also been ob¬ 
served by the OGLE collaboration (specifically MOA-BLG-2015-020 and MOA-BLG- 
2015-245 are also known as OGLE-BLG-2015-0102 and OGLE-BLG-2015-1286, respec¬ 
tively) and 64 out of the 160 OGLE events have also been observed by the MOA 
collaboration (for the detail of this we refer to the MOA microlensing alert webpage 
http://www.phys.canterbury.ac.nz/moa/microlensing_alerts.html). In Fig. 6 we show the 
corresponding lightcurves for the 3 variants of our photometry reduction for all the events 
(online material). The first method, which relies on an analysis based on the OGLE on¬ 
line catalog of nearby stars to the microlensed source, is not presented for MOA-only 
events and for 4 OGLE events (OGLE-BLG-2015-0561, OGLE-BLG-2015-0941, OGLE- 
BLG-2015-1112 and OGLE-BLG-2015-1234). 

^This event is also known as OGLE-BLG-2015-1253. 



